A Software-Defined Radio 
for the Masses, Part I 


This serves describes a complete PC-based, software-defined 
radio that uses a sound card and an innovative detector 
circuit. Mathematics 1s minimized in the 
explanation. Come see how it’s done. 


certain convergence occurs 

when multiple technologies 

align in time to make possible 
those things that once were only 
dreamed. The explosive growth of the 
Internet starting in 1994 was one of 
those events. While the Internet had 
existed for many years in government 
and education prior to that, its popu- 
larity had never crossed over into the 
general populace because of its slow 
speed and arcane interface. The devel- 
opment of the Web browser, the 
rapidly accelerating power and avail- 
ability of the PC, and the availability 
of inexpensive and _ increasingly 
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By Gerald Youngblood, AC5OG 


speedy modems brought about the 
Internet convergence. Suddenly, it all 
came together so that the Internet and 
the worldwide Web joined the every- 
day lexicon of our society. 

A similar convergence is occurring 
in radiocommunications through digi- 
tal signal processing (DSP) software to 
perform most radio functions at per- 
formance levels previously considered 
unattainable. DSP has now been 
incorporated into much of the ama- 
teur radio gear on the market to de- 
liver improved noise-reduction and 
digital-filtering performance. More 
recently, there has been a lot of discus- 
sion about the emergence of so-called 
software-defined radios (SDRs). 

A software-defined radio is charac- 
terized by its flexibility: Simply modi- 
fying or replacing software programs 


can completely change its functional- 
ity. This allows easy upgrade to new 
modes and improved performance 
without the need to replace hardware. 
SDRs can also be easily modified to 
accommodate the operating needs of 
individual applications. There is a dis- 
tinct difference between a radio that 
internally uses software for some of its 
functions and a radio that can be com- 
pletely redefined in the field through 
modification of software. The latter is 
a software-defined radio. 

This SDR convergence is occurring 
because of advances in software and 
silicon that allow digital processing of 
radio-frequency signals. Many of 
these designs incorporate mathemati- 
cal functions into hardware to perform 
allofthe digitization, frequency selec- 
tion, and down-conversion to base- 
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band. Such systems can be quite com- 
plex and somewhat out of reach to 
most amateurs. 

One problem has been that unless 
you are a math wizard and proficient 
in programming C++ or assembly lan- 
guage, you are out of luck. Each can be 
somewhat daunting to the amateur as 
well as to many professionals. Two 
years ago, I set out to attack this chal- 
lenge armed with a fascination for 
technology and a 25-year-old, virtu- 
ally unused electrical engineering de- 
gree. I had studied most of the math in 
college and even some of the signal 
processing theory, but 25 years is a 
long time. I found that it really was a 
challenge to learn many of the disci- 
plines required because much of the 
literature was written from a math- 
ematician’s perspective. 

Now that I am beginning to grasp 
many of the concepts involved in soft- 
ware radios, I want to share with the 
Amateur Radio community what I 
have learned without using much 
more than simple mathematical con- 
cepts. Further, a software radio 
should have as little hardware as pos- 
sible. If you have a PC with a sound 
card, you already have most of the 
required hardware. With as few as 
three integrated circuits youcan be up 
and running with a Tayloe detector— 
an innovative, yet simple, direct-con- 
version receiver. With less than a 
dozen chips, you can build a trans- 
ceiver that will outperform much of 
the commercial gear on the market. 


Approach the Theory 


In this article series, Ihave chosen to 
focus on practical implementation 
rather than on detailed theory. There 
are basic facts that must be understood 
to build a software radio. However, 
much like working with integrated cir- 
cuits, you don’t have to know how to 
create the IC in order to use it in a de- 
sign. The convention I have chosen is to 
describe practical applications fol- 
lowed by references where appropriate 
for more detailed study. One of the 
easier to comprehend references I have 
found is The Scientist and Engineer’s 
Guide to Digital Signal Processing by 
Steven W. Smith. Itis free for download 
over the Internet at www.DSPGuide. 
com. I consider it required reading for 
those who want to dig deeper into 
implementation aswell as theory. I will 
refer to it as the “DSP Guide” many 
times in this article series for further 
study. 

So get out your four-function calcu- 
lator (okay, maybe you need six or 
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seven functions) and let’s get started. 

But first, let’s set forth the objectives 

of the complete SDR design: 

* Keep the math simple 

¢ Use a sound-card equipped PC to pro- 
vide all signal-processing functions 

¢ Program the user interface and all 
signal-processing algorithms in 
Visual Basic for easy development 
and maintenance : 

¢ Utilize the Intel Signal Processing 
Library for core DSP routines to 
minimize the technical knowledge 
requirement and development time, 
and to maximize performance 

«Integrate a direct conversion (D-C) 
receiver for hardware design sim- 
plicity and wide dynamic range 

«Incorporate direct digital synthesis 
(DDS) to allow flexible frequency 
control 

«Include transmit capabilities using 
similar techniques as those used in 
the D-C receiver. 


Analog and Digital Signals in 
the Time Domain 


To understand DSP we first need to 
understand the relationship between 
digital signals and their analog coun- 
terparts. If we look at a 1-V (pk) sine 
wave on an analog oscilloscope, we see 
that the signal makes a perfectly 
smooth curve on the scope, no matter 
how fast the sweep frequency. In fact, 
ifit were possible to build a scope with 
an infinitely fast horizontal sweep, it 
would still display a perfectly smooth 
curve (really a straight line at that 
point). As such, it is often called acon- 
tinuous-time signal since it is continu- 
ous in time. In other words, there are 
an infinite number of different volt- 
ages along the curve, as can be seen on 
the analog oscilloscope trace. 

On the other hand, if we were to 
measure the same sine wave with a 
digital voltmeter at a sampling rate of 
four times the frequency of the sine 
wave, starting at time equals zero, we 
would read: 0 V at 0°, 1 Vat 90°, 0 V at 
180° and —1 V at 270° over one com- 
plete cycle. The signal could continue 
perpetually, and we would still read 
those same four voltages over and 
again, forever. We have measured the 
voltage of the signal at discrete mo- 
ments in time. The resulting voltage- 
measurement sequence is therefore 
called a discrete-time signal. 

If we save each discrete-time signal 
voltage in a computer memory and we 
know the frequency at which we 
sampled the signal, we havea discrete- 
time sampled signal. This is what an 
analog-to-digital converter (ADC) 


does. It uses a sampling clock to mea- 
sure discrete samples of an incoming 
analog signal at precise times, and it 
produces a digital representation of 
the input sample voltage. 

In 1933, Harry Nyquist discovered 
that to accurately recover all the com- 
ponents of a periodic waveform, it is 
necessary to use a sampling frequency 
of at least twice the bandwidth of the 
signal being measured. That mini- 
mum sampling frequency is called the 
Nyquist criterion. This may be ex- 
pressed as: 


fs 22 fow 


where f, is the sampling rate and /py, is 
the bandwidth. See? The math isn’t so 
bad, is it? 

Now as an example of the Nyquist 
criterion, let’s consider human hear- 
ing, which typically ranges from 20 Hz 
to 20 kHz. To recreate this frequency 
response, a CD player must sample at 
a frequency of at least 40 kHz. As we 
will soon learn, the maximum fre- 
quency component must be limited to 
20 kHz through low-pass filtering to 
prevent distortion caused by false im- 
ages of the signal. To ease filter re- 
quirements, therefore, CD players use 
a standard sampling rate of 44,100 Hz. 
All modern PC sound cards support 
that sampling rate. 

What happens if the sampled band- 
width is greater than halfthe sampling 
rate and is not limited by a low-pass 
filter? An alias of the signal is produced 
that appears in the output along with 
the original signal. Aliases can cause 
distortion, beat notes and unwanted 
spurious images. Fortunately, alias 
frequencies can be precisely predicted 
and prevented with proper low-pass or 
band-pass filters, which are often re- 
ferred to as anti-aliasing filters, as 
shown in Fig 1. There are even cases 
where the alias frequency can be used 
to advantage; that will be discussed 
later in the article. 

This is the point where most texts 
on DSP gointo great detail about what 
sampled signals look like above the 
Nyquist frequency. Since the goal of 
this article is practical implementa- 
tion, I refer you to Chapter 3 of the 
DSP Guide for a more in-depth discus- 
sion of sampling, aliases, A-to-D and 


(Eq 1) 


Signal 


Fig 1—A/D conversion with antialiasing 
low-pass filter. 
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D-to-A conversion. Also refer to Doug Smith’s article, “Sig- 
nals, Samples, and Stuff: A DSP Tutorial.”} 

What you need to know for now is that if we adhere to the 
Nyquist criterion in Eq 1, we can accurately sample, pro- 
cess and recreate virtually any desired waveform. The 
sampled signal will consist of a series of numbers in com- 
puter memory measured at time intervals equal to the 
sampling rate. Since we now know the amplitude of the 
signal at discrete time intervals, we can process the digi- 
tized signal in software with a precision and flexibility not 
possible with analog circuits. 


From RF to a PC’s Sound Card 


Our objective is to convert a modulated radio-frequency 
signal from the frequency domain to the time domain for 
software processing. In the frequency domain, we measure 
amplitude versus frequency (as with a spectrum analyzer); 
in the time domain, we measure amplitude versus time (as 
with an oscilloscope). 

Inthis application, we choose to use a standard 16-bit PC 
sound card that has a maximum sampling rate of 
44,100 Hz. According to Eq 1, this means that the maxi- 
mum-bandwidth signal we can accommodate is 22,050 Hz. 
With quadrature sampling, discussed later, this can actu- 
ally be extended to 44 kHz. Most sound cards have built-in 
antialiasing filters that cut off sharply at around 20 kHz. 
(For a couple hundred dollars more, PC sound cards are 
now available that support 24 bits at a 96-kHz sampling 
rate with up to 105 dB of dynamic range.) 

Most commercial and amateur DSP designs use dedicated 
DSPs that sample intermediate frequencies (IFs) of 40 kHz 
or above. They use traditional analog superheterodyne tech- 
niques for down-conversion and filtering. With the advent of 
very-high-speed and wide-bandwidth ADCs, it is now pos- 
sible to directly sample signals up through the entire HF 
range and even into the low VHF range. For example, the 
Analog Devices AD9430 A/D converter is specified with 
sample rates up to 210 Msps at 12 bits of resolution and a 
700-MHz bandwidth. That 700-MHz bandwidth can be used 
in under-sampling applications, a topic that is beyond the 
scope of this article series. 

The goal of my project is to build a PC-based software- 
defined radio that uses as little external hardware as pos- 
sible while maximizing dynamic range and flexibility. To 
do so, we will need to convert the RF signal to audio fre- 
quencies in a way that allows removal of the unwanted 
mixing products or images caused by the down-conversion 
process. The simplest way to accomplish this while main- 
taining wide dynamic range is to use D-C techniques to 
translate the modulated RF signal directly to baseband. 


‘Notes appear on page 21. 


LPF 
1.5 kHz 
14.001 MHz ——> Og Baseband 


14.0 MHz 


Fig 2—A direct-conversion real mixer with a 1.5-kHz low-pass 
filter. 
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We can mix the signal with an oscillator tuned to the RF 
carrier frequency to translate the bandwidth-limited sig- 
nal to a 0-Hz IF as shown in Fig 2. 

The example in the figure shows a 14.001-MHz carrier 
signal mixed with a 14.000-MHz local oscillator to translate 
the carrier to 1 kHz. If the low-pass filter had a cutoff of 
1.5 kHz, any signal between 14.000 MHz and 14.0015 MHz 
would be within the passband of the direct-conversion re- 
ceiver. The problem with this simple approach is that we 
would also simultaneously receive all signals between 
13.9985 MHz and 14.000 MHz as unwanted images within 
the passband, as illustrated in Fig 3. Why is that? 

Most amateurs are familiar with the concept of sum and 
difference frequencies that result from mixing two signals. 
When a carrier frequency, f,, is mixed with a local oscilla- 
tor, f\,, they combine in the general form: 


1 
fefo =5\ Fe + fio) + Fe - fio) | Ba?) 
When we use the direct-conversion mixer shown in Fig 2, 
we will receive these primary output signals: 


fe + fig =14.001 MHz +14.000 MHz = 28.001 MHz 
fe - fio = 14.001 MHz - 14.000 MHz = 0.001 MHz 


Note that we also receive the image frequency that “folds 
over” the primary output signals: 


~ fo + fig = 14.001 MHz + 14.000 MHz = -0.001 MHz 


A low-pass filter easily removes the 28.001-MHz sum 
frequency, but the —0.001-MHz difference-frequency image 
will remain in the output. This unwanted image is the 
lower sideband with respect to the 14.000-MHz carrier fre- 
quency. This would not be a problem if there were no sig- 
nals below 14.000 MHz to interfere. As previously stated, 
all undesired signals between 13.9985 and 14.000 MHz will 
translate into the passband along with the desired signals 
above 14.000 MHz. The image also results in increased 
noise in the output. 

So how can we remove the image-frequency signals? It 
can be accomplished through quadrature mixing. Phasing 
or quadrature transmitters and receivers—also called 
Weaver-method or image-rejection mixers—have existed 
since the early days of single sideband. In fact, my first 
SSB transmitter was a used Central Electronics 20A ex- 
citer that incorporated a phasing design. Phasing systems 
lost favor in the early 1960s with the advent of relatively 
inexpensive, high-performance filters. 

To achieve good opposite-sideband or image suppression, 
phasing systems require a precise balance of amplitude and 
phase between two samples of the signal that are 90° out 


Image Signals Real Signals 


-28.001 MHz -1 kHz 


1 kHz 


28.001 MHz 
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Fig 3—~Output spectrum of a real mixer illustrating the sum, 
difference and image frequencies. 
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of phase or in quadrature with each 
other—“orthogonal” is the term used 
in some texts. Until the advent of digi- 
tal signal processing, it was difficult 
to realize the level of image rejection 
performance required of modern radio 
systems in phasing designs. Since 
digital signal processing allows pre- 
cise numerical control of phase and 
amplitude, quadrature modulation 
and demodulation are the preferred 
methods. Such signals in quadrature 
allow virtually any modulation 
method to be implemented in software 
using DSP techniques. 


Give Me I and Q and I Can 
Demodulate Anything 


First, consider the direct-conversion 
mixer shown in Fig 2. When the RF sig- 
nal is converted to baseband audio us- 
ing a single channel, we can visualize 
the output as varying in amplitude 
along a single axis as illustrated in 
Fig 4. We will refer to this as the in- 
phase or I signal. Notice that its magni- 
tude varies from a positive value to a 
negative value at the frequency of the 
modulating signal. If we use a diode to 
rectify the signal, we would have cre- 
ated a simple envelope or AM detector. 

Remember that in AM envelope de- 
tection, both modulation sidebands 
carry information energy and both are 
desired at the output. Only amplitude 
information is required to fully de- 
modulate the original signal. The 
problem is that most other modulation 
techniques require that the phase of 
the signal be known. This is where 
quadrature detection comes in. If we 
delay a copy of the RF carrier by 90° to 
form a quadrature (Q) signal, we can 
then use it in conjunction with the 
original in-phase signal and the math 
we learned in middle school to deter- 
mine the instantaneous phase and 
amplitude of the original signal. 

Fig 5 illustrates an RF carrier with 
the level of the J signal plotted on the 
x-axis and that of the Q signal plotted 
on the y-axis of a plane. This is often 
referred to in the literature as a 
phasor diagram in the complex plane. 
We are now able to extrapolate the two 
signals to draw an arrow or phasor 
that represents the instantaneous 
magnitude and phase of the original 
signal. 

Okay, here is where you will have to 
use a couple of those extra functions 
on the calculator. To compute the 
magnitude m;, or envelope of the sig- 
nal, we use the geometry of right tri- 
angles. In a right triangle, the square 
of the hypotenuse is equal to the sum 
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of the squares of the other two sides— 
according to the Pythagorean theo- 
rem. Or restating, the hypotenuse as 
m, (magnitude with respect to time): 


2 2 
m =y1; +0; 


The instantaneous phase of the sig- 
nal as measured counterclockwise 
from the positive J axis and may be 
computed by the inverse tangent (or 
arctangent) as follows: 


= tannt| 2 
%, = tan 24) 


Therefore, if we measured the in- 
stantaneous values of J and Q, we 
would know everything we needed to 
know about the signal at a given mo- 
ment in time. This is true whether we 
are dealing with continuous analog 
signals or discrete sampled signals. 
With J and Q, we can demodulate AM 
signals directly using Eq 3 and FM 
signals using Eq 4. To demodulate 
SSB takes one more step. Quadrature 
signals can be used analytically to re- 
move the image frequencies and leave 
only the desired sideband. 

The mathematical equations for 
quadrature signals are difficult but 
are very understandable with a little 
study.2 I highly recommend that you 
read the online article, “Quadrature 


(Eq 3) 


(Eq 4) 


—_———— ene — [n Phase 
mt 


Fig 4—An in-phase signal (/) on the real 
plane. The magnitude, my, is easily 
measured as the instantaneous peak 
voltage, but no phase information is 
available from in-phase detection. This is 
the way an AM envelope detector works. 


Signals: Complex, But Not Compli- 
cated,” by Richard Lyons. It can be 
found at www.dspguru.com/info/ 
tutor/quadsig.htm. The article de- 
velops in a very logical manner how 
quadrature-sampling J/Q demodula- 
tion is accomplished. A basic under- 
standing of these concepts is essential 
to designing software-defined radios. 
We can take advantage of the ana- 
lytic capabilities of quadrature signals 
through a quadrature mixer. To under- 
stand the basic concepts of quadrature 
mixing, refer to Fig 6, which illustrates 
a quadrature-sampling //Q mixer. 
First, the RF input signal is band- 
pass filtered and applied to the two 
parallel mixer channels. By delaying 
the local oscillator wave by 90°, we can 
generate acosine wave that, in tandem, 
forms a quadrature oscillator. The RF 
carrier, f,(¢), is mixed with the respec- 
tive cosine and sine wave local oscilla- 
tors and is subsequently low-pass 
filtered to create the in-phase, /(t), and 
quadrature, Q(t), signals. The Q(t) 


mM 


2 


Fig 5—1 +/Q are shown on the complex 
plane. The vector rotates counterclock- 
wise at a rate of 2xf,_. The magnitude and 
phase of the rotating vector at any instant 
in ane may be determined through Eqs 3 
and 4. 


\(t) 


Q(t) 


Fig 6—Quadrature sampling mixer: The RF carrier, f,, is fed to parallel mixers. The local 
oscillator (Sine) is fed to the lower-channel mixer directly and Is delayed by 90° (Cosine) 
to feed the upper-channel mixer. The low-pass filters provide antialias filtering before 
analog-to-digital conversion. The upper channel provides the in-phase (/,)) signal and the 


lower channel provides the quadrature (Q,,) 
and A/D converters are integrated on the p 


signal. In the PC SDR the low-pass filters 
sound card. 
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channel is phase-shifted 90° relative to 
the [(t) channel through mixing with 
the sine local oscillator. The low-pass 
filter is designed for cutoff below the 
Nyquist frequency to prevent aliasing 
in the A/D step. The A/D converts con- 
tinuous-time signals to discrete-time 
sampled signals. Now that we have the 
I and Q samples in memory, we can 
perform the magic of digital signal pro- 
cessing. 

Before we go further, let me reiter- 
ate that one of the problems with this 
method of down-conversion is that it 
can be costly to get good opposite-side- 
band suppression with analog circuits. 
Any variance in component values will 
cause phase or amplitude imbalance 
between two channels, resulting in a 
corresponding decrease in opposite- 
sideband suppression. With analog 
circuits, it is difficult to achieve better 
than 40 dB of suppression without 
much higher cost. Fortunately, it is 
straightforward to correct the analog 
imbalances in software. 

Another significant drawback of di- 
rect-conversion receivers is that the 
noise increases as the demodulated sig- 
nal approaches 0 Hz. Noise contribu- 
tions come from a number of sources, 
such as 1/fnoise from the semiconduc- 
tor devices themselves, 60-Hz and 
120-Hz line noise or hum, microphonic 
mechanical noise and local-oscillator 
phase noise near the carrier frequency. 
This can limit sensitivity since most 
people prefer their CW tones to be be- 
low 1 kHz. It turns out that most of 
the low-frequency noise rolls off above 
1 kHz. Since a sound card can process 
signals all the way up to 20 kHz, why 
not use some of that bandwidth to move 
away from the low frequency noise? The 
PC SDR uses an 11.025-kHz, offset- 
baseband IF to reduce the noise to a 
manageable level. By offsetting the 
local oscillator by 11.025 kHz, we can 
now receive signals near the carrier 


Frequency 
Fig 7—FFT output resembles a comb filter: 
Each bin of the FFT overlaps its adjacent 
bins just as in a comb filter. The 3-dB 
points overlap to provide linear output. The 
phase and magnitude of the signal in each 


bin is easily determined mathematically 
with Eqs 3 and 4. 
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frequency without any of the low- 
frequency noise issues. This also 
significantly reduces the effects of lo- 
cal-oscillator phase noise. Once we 
have digitally captured the signal, it is 
a trivial software task to shift the de- 
modulated signal down to a 0-Hz offset. 


DSP in the Frequency Domain 


Every DSP text I have read thus far 
concentrates on time-domain filtering 
and demodulation of SSB signals us- 
ing finite-impulse-response (FIR) fil- 
ters. Since these techniques have been 
thoroughly discussed in the litera- 
ture!.3: 4 and are not currently used in 
my PC SDR, they will not be covered 
in this article series. 

My PC SDR uses the power of the 
fast Fourier transform (FFT) to do al- 
most all of the heavy lifting in the fre- 
quency domain. Most DSP texts use a 
lot of ink to derive the math so that one 
can write the FFT code. Since Intel has 
so helpfully provided the code in ex- 
ecutable form in their signal-process- 
ing library,5 we don’t care how to write 
an FFT: We just need to know how to 
use it. Simply put, the FFT converts 
the complex / and Q discrete-time sig- 
nals into the frequency domain. The 
FFT output can be thought of as a 
large bank of very narrow band-pass 
filters, called bins, each one measur- 
ing the spectral energy within its 
respective bandwidth. The output re- 
sembles acomb filter wherein each bin 
slightly overlaps its adjacent bins 
forming a scalloped curve, as shown in 
Fig 7. When asignalis precisely at the 
center frequency ofa bin, there will be 
a corresponding value only in that bin. 
As the frequency is offset from the 
bin’s center, there will be a corre- 
sponding increase in the value of the 


adjacent bin and a decrease in the 
value of the current bin. Mathemati- 
cal analysis fully describes the rela- 
tionship between FFT bins,® but such 
is beyond the scope of this article. 

Further, the FFT allows us to mea- 
sure both phase and amplitude of the 
signal within each bin using Eqs 3 and 
4 above. The complex version allows us 
to measure positive and negative fre- 
quencies separately. Fig 8 illustrates 
the output of a complex, or quadra- 
ture, FFT. 

The bandwidth of each FFT bin may 
be computed as shown in Eq 5, where 
BWpin is the bandwidth of a single bin, 
f,is the sampling rate and Nis the size 
of the FFT. The center frequency of 
each FFT bin may be determined by 
Eq 6 where frenter is the bin’s center 
frequency, nis the bin number, f, is the 
sampling rate and N is the size of the 
FFT. Bins zero through (N/2)-1 repre- 
sent upper-sideband frequencies and 
bins N/2 to N-1 represent lower-side- 
band frequencies around the carrier 
frequency. 


BWein = 4 (Eq 5) 
Senter = whe (Eq 6) 


If we assume the sampling rate of 
the sound card is 44.1 kHz and the 
number of FFT bins is 4096, then the 
bandwidth and center frequency of 
each bin would be: 


_ 44100 


bin ~ 4096 


= 10.7666 Hz and 


Senter = 210.7666 Hz 


What this all means is that the 
receiver will have 4096, ~11-Hz-wide 


Center = nfs 
} LSB f ino) USB } 
N N-4 ad 
2 (Bin 4095) at 
(Bin 2048) (Bin 2047) 


Fig 8—Compiex FFT output: The output of a complex FFT may be thought of as a series 
of band-pass filters aligned around the carrier frequency, f,, at bin 0. N represents the 
number of FFT bins. The upper sideband is located in bins 1 through (N/2}-1 and the 
lower sideband is located in bins N/2 to N-1. The center frequency and bandwidth of 
each bin may be calculated using Eqs 5 and 6. 
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band-pass filters. We can therefore 
create band-pass filters from 11 Hz to 
approximately 40 kHz in 11-Hz steps. 
The PC SDR performs the following 
functions in the frequency domain af- 
ter FFT conversion: 
* Brick-wall fixed and variable band- 
pass filters 
* Frequency conversion 
«SSB/CW demodulation 
¢ Sideband selection 
¢ Frequency-domain noise subtraction 
* Frequency-selective squelch 
* Noise blanking 
* Graphic equalization (“tone control”) 
«Phase and amplitude balancing to 
remove images 
«SSB generation 
¢ Future digital modes such as PSK31 
and RTTY 
Once the desired frequency-domain 
processing is completed, it is simple to 
convert the signal back to the time do- 
main by using an inverse FFT. In the 
PC SDR, only AGC and adaptive noise 
filtering are currently performed in the 
time domain. A simplified diagram of 
the PC SDR software architecture is 
provided in Fig 9. These concepts 
will be discussed in detail in a future 
article. 


Sampling RF Signals with the 
Tayloe Detector: A New Twist 
on an Old Problem 


While searching the Internet for 
information on quadrature mixing, I 
ran across a most innovative and el- 
egant design by Dan Tayloe, N7VE. 
Dan, who works for Motorola, has de- 
veloped and patented (US Patent 
#6,230,000) what has been called the 
Tayloe detector’. The beauty of the 
Tayloe detector is found in both its 
design elegance and its exceptional 
performance. It resembles other con- 
cepts in design, but appears unique in 
its high performance with minimal 
components.®: 9 10, 1! In its simplest 
form, you can build a complete quadra- 
ture down converter with only three or 
four ICs (less the local oscillator) at a 
cost of less than $10. 

Fig 10 illustrates a single-balanced 
version of the Tayloe detector. It can be 
visualized as a four-position rotary 
switch revolving at a rate equal to the 
carrier frequency. The 50-Q antenna 
impedance is connected to the rotor and 
each of the four switch positions is con- 
nected to a sampling capacitor. Since 
the switch rotor is turning at exactly 
the RF carrier frequency, each capaci- 
tor will track the carrier’s amplitude 
for exactly one-quarter of the cycle and 
will hold its value for the remainder of 
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Fig 9—SDR receiver software architecture: The / and Q signals are fed from the sound- 
card input directly to a 4096-bin complex FFT. Band-pass filter coefficients are 
precomputed and converted to the frequency domain using another FFT. The frequency- 
domain filter is then multiplied by the frequency-domain signal to provide brick-wall 
filtering. The filtered signal is then converted to the time domain using the inverse FFT. 
Adaptive noise and notch filtering and digital AGC follow in the time domain. 


Fig 10—Tayloe detector: The switch rotates at the carrier frequency so that each 
capacitor samples the signal once each revolution. The 0° and 180° capacitors 
differentially sum to provide the in-phase (/ signal and the 90° and 270° capacitors sum 
to provide the quadrature (Q) signal. 


the cycle. The rotating switch will 
therefore sample the signal at 0°, 90°, 
180° and 270°, respectively. 

As shown in Fig 11, the 50-Q imped- 
ance of the antenna and the sampling 
capacitors form an R-C low-pass filter 
during the period when each respec- 
tive switch is turned on. Therefore, 
each sample represents the integral or 
average voltage of the signal during its 
respective one-quarter cycle. When 
the switch is off, each sampling capaci- 
tor will hold its value until the next 
revolution. If the RF carrier and the 
rotating frequency were exactly in 
phase, the output of each capacitor 
will be a dc level equal to the average 


Rant 


! feat Baseband 
tee 


Fig 11—Track and hold sampling circuit: 
Each of the four sampling capacitors in the 
Tayloe detector form an RC track-and-hold 
circuit. When the switch is on, the 
capacitor will charge to the average value 
of the carrier during its respective one- 
quarter cycle. During the remaining three- 
quarters cycle, it will hold its charge. The 
local-oscillator frequency is equal to the 
carrier frequency so that the output will be 
at baseband. 
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value of the sample. 

If we differentially sum outputs of 
the 0° and 180° sampling capacitors 
with an op amp (see Fig 10), the out- 
put would be a de voltage equal to two 
times the value of the individually 
sampled values when the switch rota- 
tion frequency equals the carrier fre- 
quency. Imagine, 6 dB of noise-free 
gain! The same would be true for the 
90° and 270° capacitors as well. The 
0°/180° summation forms the J chan- 
nel and the 90°/270° summation forms 
the Q channel of the quadrature down- 
conversion. 

As we shift the frequency of the car- 
rier away from the sampling fre- 
quency, the values of the inverting 
phases will no longer be dc levels. The 
output frequency will vary according 
to the “beat” or difference frequency 
between the carrier and the switch-ro- 
tation frequency to provide an accu- 
rate representation of all the signal 


components converted to baseband. 

Fig 12 provides the schematic for a 
simple, single-balanced Tayloe detec- 
tor. It consists of a PI5V331, 1:4 FET 
demultiplexer that switches the signal 
to each of the four sampling capaci- 
tors. The 74AC74 dual flip-flop is con- 
nected as a divide-by-four Johnson 
counter to provide the two-phase clock 
to the demultiplexer chip. The outputs 
of the sampling capacitors are differ- 
entially summed through the two 
LT1115 ultra-low-noise op amps to 
form the] and Q outputs, respectively. 
Note that the impedance of the 
antenna forms the input resistance for 
the op-amp gain as shown in Eq7. This 
impedance may vary significantly 
with the actual antenna. I use instru- 
mentation amplifiers in my final de- 
sign to eliminate gain variance with 
antenna impedance. More informa- 
tion on the hardware design will be 
provided in a future article. 


Since the duty cycle of each switch 
is 25%, the effective resistance in the 
RC network is the antenna impedance 
multiplied by four in the op-amp gain 
formula, as shown in Eq 7: 


_ _ Re 
~ 4R 


(Eq 7) 


ant 


For example, with a feedback resis- 
tance, R,, of 3.3 kQ and antenna im- 
pedance, Ray, of 50 Q, the resulting 
gain of the input stage is: 


3300 
4x50 


G= =16.5 


The Tayloe detector may also be 
analyzed as a digital commutating fil- 
ter,.12,13,14 This means that it operates 
as avery-high-Q tracking filter, where 
Eq 8 determines the bandwidth and n 
is the number of sampling capacitors, 
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Fig 12—Singly balanced Tayloe detector. 
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Rant is the antenna impedance and C, 
is the value of the individual sampling 
capacitors. Eq 9 determines the Qget 
of the filter, where f, is the center fre- 
quency and BW,,; is the bandwidth of 
the filter. 


oe eee (Eq 8) 


-—Z (Eq 9) 
oat BWaet 

By example, if we assume the sam- 
pling capacitor to be 0.27 pF and the 
antenna impedance to be 50 Q, then 
BW and Q are computed as follows: 


1 


BWee, => = 5895 Hz 
(2)(4)(50)(2.7x 1077) 
6 
_, = 14.001x10 _ 9475 
5895 


Since the PC SDR uses an offset 
baseband IF, I have chosen to design 
the detector’s bandwidth to be 40 kHz 
to allow low-frequency noise elimina- 
tion as discussed above. 

The real payoff in the Tayloe detec- 
tor is its performance. It has been 
stated that the ideal commutating 
mixer has a minimum conversion loss 
(which equates to noise figure) of 
3.9 dB.15. 16 Typical high-level diode 
mixers have aconversion loss of 6-7 dB 
and noise figures 1 dB higher than the 
loss. The Tayloe detector has less than 
1 dB of conversion loss, remarkably. 
How can this be? The reason is that it 
is not really a mixer but a sampling 
detector in the form of a quadrature 
track and hold. This means that the 
design adheres to discrete-time sam- 
pling theory, which, while similar to 
mixing, has its own unique character- 
istics. Because a track and hold actu- 
ally holds the signal value between 
samples, the signal output never goes 
to zero. 

This is where aliasing can actually 
be used to our benefit. Since each 
switch and capacitor in the Tayloe 
detector actually samples the RF sig- 
nal once each cycle, it will respond to 
alias frequencies as well as those 
within the Nyquist frequency range. 
In a traditional direct-conversion re- 
ceiver, the local-oscillator frequency is 
set to the carrier frequency so that the 
difference frequency, or IF, is at 0 Hz 
and the sum frequency is at two times 
the carrier frequency per Eq 2. We 
normally remove the sum frequency 
through low-pass filtering, resulting 
in conversion loss and a corresponding 
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increase in noise figure. In the Tayloe 
detector, the sum frequency resides at 
the first alias frequency as shown in 
Fig 13. Remember that an alias is a 
real signal and will appear in the out- 
put as if it were a baseband signal. 
Therefore, the alias adds to the base- 
band signal for a theoretically loss- 
less detector. In real life, there is a 
slight loss due to the resistance of the 
switch and aperture loss due to imper- 
fect switching times. 


PC SDR Transceiver Hardware 


The Tayloe detector therefore pro- 
vides a low-cost, high-performance 
method for both quadrature down-con- 
version as well as up-conversion for 
transmitting. For a complete system, 
we would need to provide analog AGC 
to prevent overload of the ADC inputs 
and a means of digital frequency con- 
trol. Fig 14 illustrates the hardware 


architecture of the PC SDR receiver as 
it currently exists. The challenge has 
been to build a low-noise analog chain 
that matches the dynamic range of the 
Tayloe detector to the dynamic range 
of the PC sound card. This will be cov- 
ered in a future article. 

I am ccurrently prototyping a 
complete PC SDR transceiver, the 
SDR-1000, that will provide general- 
coverage receive from 100 kHz to 
54 MHz and will transmit on all ham 
bands from 160 through 6 meters. 


SDR Applications 


At the time of this writing, the typi- 
cal entry-level PC now runs at a clock 
frequency greater than 1 GHz and 
costs only a few hundred dollars. We 
now have exceptional processing 
power at our disposal to perform DSP 
tasks that were once only dreams. The 
transfer of knowledge from the aca- 


Fig 13—Alias summing on Tayloe detector output: Since the Tayloe detector samples the 
signal the sum frequency (fc + fs) and its image (-fc — fs) are located at the first alias 
frequency. The alias signals sum with the baseband signals to eliminate the mixing 
product loss associated with traditional mixers. In a typical mixer, the sum frequency 
energy is lost through filtering thereby increasing the noise figure of the device. 
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Fig 14—PC SDR receiver hardware architecture: After band-pass filtering the antenna is 
fed directly to the Tayloe detector, which in turn provides {and Q outputs at baseband. A 
DDS and a divide-by-four Johnson counter drive the Tayloe detector demultiplexer. The 
LT1115s offer ultra-low noise-differential summing and amplification prior to the wide- 
dynamic-range analog AGC circuit formed by the SSM2164 and AD8307 log amplifier. 
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demic to the practical is the primary 
limit of the availability of this technol- 
ogy to the Amateur Radio experi- 
menter. This article series attempts to 
demystify some of the fundamental 
concepts to encourage experimenta- 
tion within our community. The ARRL 
recently formed aSDR Working Group 
for supporting this effort, as well. 
The SDR mimics the analog world in 
digital data, which can be manipu- 
lated much more precisely. Analog 
radio has always been modeled math- 
ematically and can therefore be pro- 
cessed in a computer. This means that 
virtually any modulation scheme may 
be handled digitally with performance 
levels difficult, or impossible, to attain 
with analog circuits. Let’s consider 
some of the amateur applications for 
the SDR: 
« Competition-grade HF transceivers 
* High-performance IF for microwave 
bands 
«Multimode digital transceiver 
* EME and weak-signal work 
* Digital-voice modes 
* Dream it and code it 


For Further Reading 


For more in-depth study of DSP 
techniques, I highly recommend that 
you purchase the following texts in 
order of their listing: 

Understanding Digital Signal Pro- 
cessing by Richard G. Lyons (see Note 
6). This is one of the best-written text- 
books about DSP. 

Digital Signal Processing Technol- 
ogy by Doug Smith (see Note 4). This 
new book explains DSP theory and 
application from an Amateur Radio 
perspective. 

Digital Signal Processing in Com- 
munications Systems by Marvin E. 
Frerking (see Note 3). This book re- 
lates DSP theory specifically to modu- 
lation and demodulation techniques 
for radio applications. 
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Where Do We Go From Here? 


Three future articles will describe 
the construction and programming of 
the PC SDR. The next article in the 
series will detail the software interface 
to the PC sound card. Integrating full- 
duplex sound with DirectX was one of 
the more challenging parts of the 
project. The third article will describe 
the Visual Basic code and the use of the 
Intel Signal Processing Library for 
implementing the key DSP algorithms 
in radio communications. The final 
article will describe the completed 
transceiver hardware for the SDR- 
1000. 
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A Software-Defined Radio 
for the Masses, Part 2 


Come learn how to use a PC sound card to enter 
the wonderful world of digital signal processing. 


art 1 gave a general description 

of digital signal processing 

(DSP) in software-defined ra- 
dios (SDRs).! It also provided an over- 
view of a full-featured radio that uses 
a personal computer to perform all 
DSP functions. This article begins de- 
sign implementation with a complete 
description of software that provides 
a full-duplex interface to a standard 
PC sound card. 

To perform the magic of digital sig- 
nal processing, we must be able to con- 
vert a signal from analog to digital and 
back to analog again. Most amateur 
experimenters already have this ca- 


1Notes appear on page 18. 
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By Gerald Youngblood, AC5OG 


pability in their shacks and many 
have used it for slow-scan television 
or the new digital modes like PSK31. 

Part 1 discussed the power of 
quadrature signal processing using in- 
phase (I) and quadrature (Q) signals 
to receive or transmit using virtually 
any modulation method. Fortunately, 
all modern PC sound cards offer the 
perfect method for digitizing the J and 
Q signals. Since virtually all cards to- 
day provide 16-bit stereo at 44-kHz 
sampling rates, we have exactly what 
we need capture and process the sig- 
nals in software. Fig 1 illustrates a 
direct quadrature-conversion mixer 
connection to a PC sound card. 

This article discusses complete 
source code for a DirectX sound-card 
interface in Microsoft Visual Basic. 
Consequently, the discussion assumes 
that the reader has some fundamen- 


tal knowledge of high-level language 
programming. 


Sound Card and PC Capabilities 


Very early PC sound cards were low- 
performance, 8-bit mono versions. To- 
day, virtually all PCs come with 
16-bit stereo cards of sufficient quality 
to be used in a software-defined radio. 
Such a card will allow us to demodu- 
late, filter and display up to approxi- 
mately a 44-kHz bandwidth, assuming 
a 44-kHz sampling rate. (The band- 
width is 44 kHz, rather than 22 kHz, 
because the use of two channels effec- 
tively doubles the sampling rate—Ed. ) 
For high-performance applications, it is 
important to select a card that offers a 
high dynamic range—on the order of 
90 dB. If you are just getting started, 
most PC sound cards will allow you to 
begin experimentation, although they 
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may offer lower performance. 

The best 16-bit price-to-perfor- 
mance ratio I have found at the time 
of this article is the Santa Cruz 6- 
channel DSP Audio Accelerator from 
Turtle Beach Inc (www.tbeach.com), 
It offers four 18-bit internal analog- 
to-digital (A/D) input channels and six 
20-bit digital-to-analog (D/A) output 
channels with sampling rates up to 
48 kHz. The manufacturer specifies a 
96-dB signal-to-noise ratio (SNR) and 
better than —-91 dB total harmonic dis- 
tortion plus noise (THD+N). Crosstalk 
is stated to be -105 dB at 100 Hz. The 
Santa Cruz card can be purchased 
from online retailers for under $70. 

Each bit on an A/D or D/A converter 
represents 6 dB of dynamic range, so 
a 16-bit converter has a theoretical 
limit of 96 dB. A very good converter 
with low-noise design is required to 
achieve this level of performance. 
Many 16-bit sound cards provide no 
more than 12-14 effective bits of dy- 
namic range. To help achieve higher 
performance, the Santa Cruz card uses 
an 18-bit A/D converter to deliver 
the 96 dB dynamic range (16-bit) 
specification. 

A SoundBlaster 64 also provides 
reasonable performance on the order 
of 76 dB SNR according to PC AV Tech 
at www.pceavtech.com. I have used 
this card with good results, but I much 
prefer the Santa Cruz card. 

The processing power needed from 
the PC depends greatly on the signal 
processing required by the application. 
Since I am using very-high-perfor- 
mance filters and large fast-Fourier 
transforms (FFTs), my applications 
require at least a 400-MHz Pentium 
II processor with a minimum of 
128 MB of RAM. If you require less 
performance from the software, you 
can get by with a much slower ma- 
chine. Since the entry level for new 
PCs is now 1 GHz, many amateurs 
have ample processing power avail- 
able. 


Microsoft DirectX versus 
Windows Multimedia 

Digital signal processing using a PC 
sound card requires that we be able to 
capture blocks of digitized J and Q data 
through the stereo inputs, process those 
signals and return them to the sound- 
card outputs in pseudo real time. This 
is called full duplex. Unfortunately, 
there is no high-level software interface 
that offers the capabilities we need for 
the SDR application. 

Microsoft now provides two appli- 
cation programming interfaces? (APIs) 
that allow direct access to the sound 
card under C++ and Visual Basic. The 
original interface is the Windows Mul- 
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timedia system using the Waveform 
Audio API. While my early work was 
done with the Waveform Audio API, I 
later abandoned it for the higher per- 
formance and simpler interface 
DirectX offers. The only limitation I 
have found with DirectX is that it does 
not currently support sound cards 
with more than 16-bits of resolution. 
For 24-bit cards, Windows Multimedia 
is required. While the Santa Cruz card 
supports 18-bits internally, it presents 
only 16-bits to the interface. For in- 
formation on where to download the 
DirectX software development kit 
(SDK) see Note 2. 


Circular Buffer Concepts 
A typical full-duplex PC sound card 
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allows the simultaneous capture and 
playback of two or more audio chan- 
nels (stereo). Unfortunately, there is 
no high-level code in Visual Basic or . 
C++ to directly support full duplex as 
required in an SDR. We will therefore 
have to write code to directly control 
the card through the DirectX API. 
DirectX internally manages all low- 
level buffers and their respective 
interfaces to the sound-card hard- 
ware. Our code will have to manage 
the high-level DirectX buffers 
(called DirectSoundBuffer and 
DirectSoundCaptureBuffer) to pro- 
vide uninterrupted operation in 
a multitasking system. The Direct- 
SoundCaptureBuffer stores the digi- 
tized signals from the stereo 
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Fig 1—Direct quadrature conversion mixer to sound-card interface used in the author's 
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A/D converter in a circular buffer and 
notifies the application upon the 
occurrence of predefined events. Once 
captured in the buffer, we can read 
the data, perform the necessary modu- 
lation or demodulation functions us- 
ing DSP and send the data to the 
DirectSoundBuffer for D/A conversion 
and output to the speakers or trans- 
mitter. 

To provide smooth operation in a 
multitracking system without audio 
popping or interruption, it will be nec- 
essary to provide a multilevel buffer for 
both capture and playback. You may 
have heard the term double buffering. 
We will use double buffering in the 
DirectSoundCaptureBuffer 
and quadruple buffering in the 
DirectSoundBuffer. I found that the 
quad buffer with overwrite detection 
was required on the output to prevent 
overwriting problems when the system 
is heavily loaded with other applica- 
tions. Figs 2A and 2B illustrate the 
concept of a circular double buffer, 
which is used for the Direct- 
SoundCaptureBuffer. Although the 
buffer is really a linear array in 
memory, as shown in Fig 2B, we can 
visualize it as circular, as illustrated in 
Fig 2A. This is so because DirectX man- 
ages the buffer so that as soon as each 
cursor reaches the end of the array, the 
driver resets the cursor to the begin- 
ning of the buffer. 

The DirectSoundCaptureBuffer is 
broken into two blocks, each equal in 
size to the amount of data to be cap- 
tured and processed between each 
event. Note that an event is much like 
an interrupt. In our case, we will use 
a block size of 2048 samples. Since we 
are using a stereo (two-channel) board 
with 16 bits per channel, we will be 
capturing 8192 bytes per block (2048 
samples x 2 channels x 2 bytes). There- 
fore, the DirectSoundCaptureBuffer 
will be twice as large (16,384 bytes). 

Since the DirectSoundCapture 
Buffer is divided into two data blocks, 
we will need to send an event notifica- 
tion to the application after each block 
has been captured. The DirectX driver 
maintains cursors that track the posi- 
tion of the capture operation at all 
times. The driver provides the means 
of setting specific locations within the 
buffer that cause an event to trigger, 
thereby telling the application to re- 
trieve the data. We may then read the 
correct block directly from the 
DirectSoundCaptureBuffer segment 
that has been completed. 

Referring again to Fig 2A, the two 
cursors resemble the hands on a clock 
face rotating in a clockwise direction. 
The capture cursor, |Play, represents 
the point at which data are currently 
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being captured. (I know that sounds 
backward, but that is how Microsoft 
defined it.) The read cursor, lWrite, 
trails the capture cursor and indicates 
the point up to which data can safely 
be read. The data after lWrite and up 
to and including 1Play are not neces- 
sarily good data because of hardware 
buffering. We can use the lWrite cur- 
sor to trigger an event that tells the 
software to read each respective block 
of data, as will be discussed later in 
the article. We will therefore receive 
two events per revolution of the circu- 
lar buffer. Data can be captured into 
one half of the buffer while data are 
being read from the other half. 

Fig 2C illustrates the Direct- 
SoundBuffer, which is used to output 
data to the D/A converters. In this case, 
we will use a quadruple buffer to allow 
plenty of room between the currently 
playing segment and the segment be- 
ing written. The play cursor, |Play, al- 
ways points to the next byte of data to 
be played. The write cursor, ]Write, is 
the point after which it is safe to write 
data into the buffer. The cursors may 
be thought of as rotating in a clockwise 
motion just as the capture cursors do. 
We must monitor the location of the 
cursors before writing to buffer loca- 
tions between the cursors to prevent 


overwriting data that have already 
been committed to the hardware for 
playback. 

Now let’s consider how the data 
maps from the DirectSoundCapture- 
Buffer to the DirectSoundBuffer. To 
prevent gaps or pops in the sound due 
to processor loading, we will want to 
fill the entire quadruple buffer before 
starting the playback looping. DirectX 
allows the application to set the start- 
ing point for the lPlay cursor and to 
start the playback at any time. 
Fig 3 shows how the data blocks map 
sequentially from the Direct- 
SoundCaptureBuffer to the Direct- 
SoundBuffer. Block 0 from the 
DirectSoundCaptureBuffer is trans- 
ferred to Block O of the Direct- 
SoundBuffer. Block 1 of the 
DirectSoundCaptureBuffer is next 
transferred to Block 1 of the 
DirectSoundBuffer and so forth. The 
subsequent source-code examples show 
how control of the buffers is accom- 
plished. 


Full Duplex, Step-by-Step 


The following sections provide a 
detailed discussion of full-duplex 
DirectX implementation. The example 
code captures and plays back a stereo 
audio signal that is delayed by four 


DirectAnimation Library 
DirectShowStream 1.0 Type Library 


Cdtcint 1.0 Type Library 
CDTCSeryv 1.0 Type Library 
[CIDXTMsft 1.0 Type Library 
CIDXTMsft3 1.0 ee Library 
1 .O.VALG Bi 


DirectSoundCaptureBuffer 


Fig 3—Method for mapping the 
DirectSoundCaptureBufter to 
the DirectSoundBufter. 


DirectSoundBuffer 


Fig ¢_Reaisiration of the DirectX8 for Visual Biber Library i in the Visual 


Basic IDE. 
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capture periods through buffering. You 
should refer to the “DirectX Audio” 
section of the DirectX 8.0 Program- 
mers Reference that is installed with 
the DirectX software developer’s kit 
(SDK) throughout this discussion. The 
DSP code will be discussed in the next 
article of this series, which will dis- 
cuss the modulation and demodula- 
tion of quadrature signals in the SDR. 
Here are the steps involved in creat- 
ing the DirectX interface: 

e Install DirectX runtime and SDK. 


e Add a reference to DirectX8 for 
Visual Basic Type Library. 

e Define Variables, I/O buffers and 
DirectX objects. 

e Implement DirectX8 events and 
event handles. 

e Create the audio devices. 

e Create the DirectX events. 

e Start and stop capture and play buff- 
ers. 

e Process the DirectXEvent8. 

e Fill the play buffer before starting 
playback. 


e Detect and correct overwrite errors. 

e Parse the stereo buffer into J and Q 
signals. 

¢ Destroy objects and events on exit. 
Complete functional source code for 

the DirectX driver written in Microsoft 

Visual Basic is provided for download 

from the QEX Web site. 


Install DirectX and Register it 
within Visual Basic 


The first step is to download the 
DirectX driver and the DirectX SDK 


Option Explicit 


‘Define Constants 

Const Fs As Long = 44100 
Const NFFT As Long = 4096 
Const BLKSIZE As Long = 2048 
Const CAPTURESIZE As Long = 


‘Define DirectX Objects 
Dim dx As New Directx8 

Dim ds As DirectSounds 

Dim 
Dim 
Dim 
Dim 


‘Define Type Definitions 
Dim dscbd As DSCBUFFERDESC 
Dim dsbd As DSBUFFERDESC 
Dim dspbd As WAVEFORMATEX 
Dim CapCurs As DSCURSORS 
Dim PlyCurs As DSCURSORS 


‘Create I/O Sound Buffers 


Dim inBuffer(CAPTURESIZE) As Integer 
Dim outBuffer(CAPTURESIZE) As Integer 


‘Define pointers and counters 
Dim Pass As Long 

Dim InPtr As Long 

Dim OutPtr As Long 

Dim StartAddr As Long 

Dim EndAddr As Long 

Dim CaptureBytes As Long 


4096 


dspb As DirectSoundPrimaryBuffer8 
dsc As DirectSoundCaptures 
dsb As DirectSoundSecondaryBuffers ‘Output Buffer object 
dscb As DirectSoundCaptureBuffer8 


‘Sampling frequency Hz 

‘Number of FFT bins 

‘Capture/play block size 
‘Capture Buffer size 


‘DirectX object 
‘DirectSound object 


‘Capture object 


‘Capture buffer description 


‘Primary buffer object 


‘Capture Buffer object 


‘DirectSound buffer description 


‘Primary buffer description 
‘DirectSound Capture Cursor 
‘DirectSound Play Cursor 


‘Number of capture passes 


‘Capture Buffer block pointer 


‘Output Buffer block pointer 


‘Buffer block starting address 


‘Ending buffer block address 
‘Capture bytes to read 


‘Demodulator Input Buffer 
‘Demodulator Output Buffer 


‘Define loop counter variables for timing the capture event cycle 
‘Start time for DirectX8Event loop 


Dim TimeStart As Double 
Dim TimeEnd As Double 
Dim AvgCtr As Long 

Dim AvgTime As Double 


‘Set up Event variables for the Capture 


Implements DirectXEvent8 
Dim hEvent(1) As Long 


Dim EVNT(1) As DSBPOSITIONNOTIFY 


Dim Receiving As Boolean 
Dim FirstPass As Boolean 


‘Ending 
‘Counts 
‘Stores 


Buffer 
DirectX Events 
for DirectX Event 


‘Allows 
‘Handle 


‘In Receive mode if true 


‘Notify position array 


time for DirectxXsEvent loop 
number of events to average 
the average event cycle time 


‘Denotes first pass from Start 


Fig 5—Declaration of variables, buffers, events and objects. This code is located in the General section of the module or form. 


88 


QEX-— Sept/Oct 2002 13 


from the Microsoft Web site (see Note 
3). Once the driver and SDK are in- 
stalled, you will need to register the 
DirectX8 for Visual Basic Type Li- 
brary within the Visual Basic devel- 
opment environment. 

If you are building the project from 
scratch, first create a Visual Basic 
project and name it “Sound.” When the 
project loads, go to the Project Menu/ 
References, which loads the form 
shown in Fig 4. Scroll through Avail- 
able References until you locate the 


DirectX8 for Visual Basic Type Library 
and check the box. When you press 
“OK,” the library is registered. 


Define Variables, Buffers and 
DirectX Objects 


Name the form in the Sound project 
frmSound. In the General section of 
frmSound, you will need to declare all 
of the variables, buffers and DirectX 
objects that will be used in the driver 
interface. Fig 5 provides the code that 
is to be copied into the General sec- 


tion. All definitions are commented in 
the code and should be self-explana- 
tory when viewed in conjunction with 
the subroutine code. 


Create the Audio Devices 


We are now ready to create the 
DirectSound objects and set up the 
format of the capture and play buff- 
ers. Refer to the source code in Fig 6 
during the following discussion. 

The first step is to create the 
DirectSound and DirectSoundCapture 


‘Set up the DirectSound Objects and the Capture and Play Buffers 


Sub CreateDevices () 
On Local Error Resume Next 


Set ds = 
Set dsc = 


dx .Direct SoundCreate (vbNullString) 
dx .DirectSoundCaptureCreate (vbNullString) 


‘Check to se if Sound Card is properly installed 


If Err.Number <> 0 Then 


‘DirectSound object 
‘DirectSound Capture 


MsgBox “Unable to start DirectSound. Check proper sound card installation” 


End 
End If 


‘Set the cooperative level to allow the Primary Buffer format to be set 
ds .SetCooperativeLevel Me.hWnd, DSSCL_PRIORITY 


‘Set up format for capture 
With dscbd 


buffer 


‘Buffer Size 


With .£xFormat 
-nFormatTag = WAVE _FORMAT_PCM 
-nChannels = 2 ‘Stereo 
.1SamplesPerSec = Fs ‘Sampling rate in Hz 
-nBitsPerSample = 16 ‘16 bit samples 
-nBlockAlign = .nBitsPerSample / 8 * .nChannels 
.lAvgBytesPerSec = .1SamplesPerSec * .nBlockAlign 
End With 
.1Flags = DSCBCAPS_DEFAULT 
.1BufferBytes = (dscbd.fxFormat.nBlockAlign * CAPTURESIZE) 
CaptureBytes = .1BufferBytes \ 2 ‘Bytes for 1/2 of capture buffer 
End With 
Set dscb = dsc.CreateCaptureBuf fer (dscbd) 


* Set up format for secondary playback buffer 


With dsbd 
.£xFormat = 
.LBufferBytes = 
.lFlags = 

End With 


aspbd = dsbd.fxFormat 
dspb.SetFormat dspbd 


Set dsb = 


End Sub 


ds .CreateSoundBuf fer (dsbd) 


dscbd.fxFormat 
dscbhd.1BufferBytes * 2 
DSBCAPS_GLOBALFOCUS Or DSBCAPS_GETCURRENTPOSITION2 


‘Create the capture buffer 


‘Play is 2X Capture Buffer Size 


‘Set Primary Buffer format 
‘to same as Secondary Buffer 


‘Create the secondary buffer 


Fig 6—Create the DirectX capture and playback devices. 
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objects. We then check for an error to 
see if we have a compatible sound card 
installed. If not, an error message would 
be displayed to the user. Next, we set 
the cooperative level DSSCL_ PRIOR- 
ITY to allow the Primary Buffer format 
to be set to the same as that of the Sec- 
ondary Buffer. The code that follows sets 
up the DirectSoundCaptureBuffer- 


‘Set events for 
Sub SetEvents() 


hEvent(0) = dx.CreateEvent (Me) 
hEvent(1) = dx.CreateEvent (Me) 
‘Buffer 


EVNT(0) .hEventNotify = 
EVNT(0).10ffset = 


‘Buffer 
EVNT(1) .hEventNotify = 
EVNT(1).10ffset = 


dscb.SetNotificationPositions 2, 


End Sub 


Description format and creates the 
DirectSoundCaptureBuffer object. The 
format is set to 16-bit stereo at the sam- 
pling rate set by the constant Fs. 
Next, the DirectSoundBuffer- 
Description is set to the same format 
as the DirectSoundCaptureBuffer- 
Description. We then set the Primary 
Buffer format to that of the Second- 


capture buffer notification at 0 and 1/2 


ary Buffer before creating the 
DirectSoundBuffer object. 


Set the DirectX Events 


As discussed earlier, the 
DirectSoundCaptureBuffer is divided 
into two blocks so that we can read 
from one block while capturing to the 
other. To do so, we must know when 


‘Event handle for first half of buffer 
‘Event handle for second half of buffer 


Event 0 sets Write at 50% of buffer 
hEvent (0) 
(dscbd.lBufferBytes \ 2) 


Event 1 Write at 100% of buffer 
hEvent (1) 
dscbhd.lBufferBytes - 1 


EVNT () 


- 1'Set event to first half of capture buffer 


‘Set Event to second half of capture buffer 


‘Set number of notification positions to 2 


a 


Fig 7—Create the DirectX events. 


a 


‘Create Devices and Set the DirectX8Events 


Private Sub Form_Load() 
CreateDevices 
SetEvents 

End Sub 


‘Create DirectSound devices 


‘Set up DirectX events 


‘Shut everything down and close application 
Private Sub Form_Unload(Cancel As Integer) 


If Receiving = True Then 


dsb.Stop 
dscb.Stop 
End If 


Dim i As Integer 


‘Stop Playback 
‘Stop Capture 


‘Kill DirectX Events 


If hEvent(i) Then dx.DestroyEvent hEvent (i) 


For i = 0 To UBound (hEvent) 
DoEvents 
Next 
Set dx = Nothing 
Set ds = Nothing 
Set dsc = Nothing 
Set dsb = Nothing 
Set dscb = Nothing 
Unload Me 
End Sub 


‘Destroy DirectX objects 


Fe A — ea ee es SSS 
Fig 8—Create and destroy the DirectSound Devices and events. 
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DirectX has finished writing to a 
block. This is accomplished using the 
DirectXEvent8. Fig 7 provides the 
code necessary to set up the two events 
that occur when the ]Write cursor has 
reached 50% and 100% of the 
DirectSoundCaptureBuffer. 

We begin by creating the two event 
handles hEvent(0) and hEvent(1). The 
code that follows creates a handle for 
each of the respective events and sets 
them to trigger after each half of the 
DirectSoundCaptureBuffer is filled. 
Finally, we set the number of notifica- 
tion positions to two and pass the 
name of the EVNT() event handle ar- 
ray to DirectX. 

The CreateDevices and SetEvents 
subroutines should be called from the 
Form_Load() subroutine. The Form_ 
Unload subroutine must stop capture 
and playback and destroy all of the 
DirectX objects before shutting down. 
The code for loading and unloading is 
shown in Fig 8. 


Starting and Stopping 
Capture/Playback 


Fig 9 illustrates how to start and 
stop the DirectSoundCaptureBuffer. 
The dscb.Start DSCBSTART_ LOOP- 
ING command starts the Direct- 
SoundCaptureBuffer in a continuous 
circular loop. When it fills the first half 
of the buffer, it triggers the DirectX 
Event8 subroutine so that the data 
can be read, processed and sent to the 
DirectSoundBuffer. Note that the 
DirectSoundBuffer has not yet been 
started since we will quadruple buffer 
the output to prevent processor load- 
ing from causing gaps in the output. 
The FirstPass flag tells the event to 
start filling the DirectSoundBuffer for 
the first time before starting the buffer 
looping. 


Processing the Direct-XEvent8 


Once we have started the Direct- 
SoundCaptureBuffer looping, the 
completion of each block will cause the 
DirectX Event8 code in Fig 10 to be 
executed. As we have noted, the events 
will occur when 50% and 100% of the 
buffer has been filled with data. Since 
the buffer is circular, it will begin 
again at the 0 location when the buffer 
is full to start the cycle all over again. 
Given a sampling rate of 44,100 Hz 
and 2048 samples per capture block, 
the block rate is calculated to be 
44,100/2048 = 21.53 blocks/s or one 
block every 46.4 ms. Since the quad 
buffer is filled before starting playback 
the total delay from input to output is 
4x 46.4 ms = 185.6 ms. 

The DirectX Event8_DXCallback 
event passes the eventid as a variable. 
The case statement at the beginning of 
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the code determines from the eventid, 
which half of the DirectSoundCapture- 
Buffer has just been filled. With that 
information, we can calculate the start- 
ing address for reading each block from 
the DirectSoundCaptureBuffer to the 
inBuffer() array with the dscb. 
ReadBuffer command. Next, we simply 
pass the inBuffer() to the external DSP 
subroutine, which returns the pro- 
cessed data in the outBuffer() array. 

Then we calculate the StartAddr 
and EndAddr for the next write loca- 
tion in the DirectSoundBuffer. Before 
writing to the buffer, we first check to 
make sure that we are not writing 
between the Write and 1Play cursors, 
which will cause portions of the buffer 
to be overwritten that have already 
been committed to the output. This 
will result in noise and distortion in 
the audio output. If an error occurs, 
the FirstPass flag is set to true and 
the pointers are reset to zero so that 
we flush the DirectSoundBuffer and 
start over. This effectively performs an 
automatic reset when the processor is 
overloaded, typically because of graph- 
ics intensive applications running 
alongside the SDR application. 

If there are no overwrite errors, we 
write the outBuffer() array that was 
returned from the DSP routine to the 
next StartAddr to EndAddr in the 
DirectSoundBuffer. Important note: In 
the sample code, the DSP subroutine 
call is commented out and the 
inBuffer() array is passed directly to 
the DirectSoundBuffer for testing of 
the code. When the FirstPass flag is 
set to True, we capture and write four 
data blocks before starting playback 
looping with the .SetCurrentPosition 
0 and .Play DSBPLAY_LOOPING 
commands. 

The subroutine calls to StartTimer 
and StopTimer allow the average com- 
putational time of the event loop to be 
displayed in the immediate window. 
This is useful in measuring the effi- 


‘Turn Capture/Playback On 
Private Sub cmdOn_Click() 
dscb.Start DSCBSTART_LOOPING 


Receiving = True 
FirstPass = True 
Start 
OutPtr = 0 
End Sub 


‘Turn Capture/Playback Off 
Private Sub cmdOff_Click() 


Receiving = False 
FirstPass = False 
dscb.Stop 
dsb.Stop 

End Sub 


ciency of the DSP subroutine code that 
is called from the event. In normal 
operation, these subroutine calls 
should be commented out. 


Parsing the Stereo Buffer 
into I and Q Signals 


One more step that is required to 
use the captured signal in the DSP 
subroutine is to separate or parse the 
left and right channel data into the J 
and @ signals, respectively. This can 
be accomplished using the code in 
Fig 11. In 16-bit stereo, the left and 
right channels are interleaved in the 
inBuffer() and outBuffer(). The code 
simply copies the alternating 16-bit 
integer values to the RealIn()), (same 
as I) and ImagIn(), (same as Q) buff- 
ers respectively. Now we are ready to 
perform the magic of digital signal 
processing that we will discuss in the 
next article of the series. 


Testing the Driver 


To test the driver, connect an audio 
generator—or any other audio device, 
such as a receiver—to the line input of 
the sound card. Be sure to mute line- 
in on the mixer control panel so that 
you will not hear the audio directly 
through the operating system. You can 
open the mixer by double clicking on 
the speaker icon in the lower right cor- 
ner of your Windows screen. It is also 
accessible through the Control Panel. 

Now run the Sound application and 
press the On button. You should hear 
the audio playing through the driver. 
It will be delayed about 185 ms from 
the incoming audio because of the qua- 
druple buffering. You can turn the 
mute contro! on the line-in mixer on 
and off to test the delay. It should 
sound like an echo. If so, you know that 
everything is operating properly. 


Coming Up Next 


In the next article, we will discuss 
in detail the DSP code that provides 


‘Start Capture Looping 
‘Set flag to receive mode 
‘This is the first pass after 


‘Starts writing to first buffer 


‘Reset Receiving flag 
‘Reset FirstPass flag 
‘Stop Capture Loop 
‘Stop Playback Loop 


Fig 9—Start and stop the capture/playback buffers. 
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‘Process the Capture events, call DSP routines, and output to Secondary Play Buffer 
Private Sub DirectXEvent8 DXCallback (ByVal eventid As Long) 


StartTimer ‘Save loop start time 
Select Case eventid ‘Determine which Capture Block is ready 
Case hEvent (0) 
InPtr = 0 ‘First half of Capture Buffer 
Case hEvent (1) 
InPtr = 1 ‘Second half of Capture Buffer 


End Select 


StartAddr = InPtr * CaptureBytes ‘Capture buffer starting address 


‘Read from DirectX circular Capture Buffer to inBuffer 
dscb.ReadBuffer StartAddr, CaptureBytes, inBuffer(0), DSCBLOCK_DEFAULT 


‘DSP Modulation/Demodulation - NOTE: THIS I§ WHERE THE DSP CODE IS CALLED 
é DSP inBuffer, outBuffer 


StartAddr = OutPtr * CaptureBytes ‘Play buffer starting address 
EndAddr = OutPtr + CaptureBytes - 1 ‘Play buffer ending address 


With dsb ‘Reference DirectSoundBuffer 
-GetCurrentPosition PlyCurs ‘Get current Play position 


‘If true the write is overlapping the lWrite cursor due to processor loading 
If PlyCurs.1lWrite >= StartAddr _ 
And PlyCurs.1Write <= EndAddr Then 
FirstPass = True ‘Restart play buffer 
OutPtr = 0 
StartAddr = 0 


End If 


‘If true the write is overlapping the 1Play cursor due to processor loading 
If PlyCurs.1Play >= StartAddr _ 
And PlyCurs.1Play <= EndAddr Then 
FirstPass = True ‘Restart play buffer 
OutPtr = 0 
StartAddr = 0 
End If 


‘Write outBuffer to DirectX circular Secondary Buffer. NOTE: writing inBuffer causes 
direct pass through. Replace 

‘with outBuffer below to when using DSP subroutine for modulation/demodulation 

-WriteBuffer StartAddr, CaptureBytes, inBuffer(0), DSBLOCK_DEFAULT 


OutPtr = IIf(OutPtr >= 3, 0, OutPtr + 1) ‘Counts 0 to 3 
If FirstPass = True Then ‘On FirstPass wait 4 counts before starting 
Pass = Pass + 1 ‘the Secondary Play buffer looping at 0 
If Pass = 3 Then ‘This puts the Play buffer three Capture cycles 
FirstPass = False ‘after the current one 
Pass = 0 ‘Reset the Pass counter 
-SetCurrentPosition 0 ‘Set playback position to zero 
.-Play DSBPLAY_ LOOPING ‘Start playback looping 
End If 
End If 
End With 
StopTimer ‘Display average loop time in immediate window 


End Sub Fig 10—Process the DirectXEvent8 event. Note that the example code passes the inBuffer() directly to the DirectSoundBuffer 
without processing. The DSP subroutine call has been commented out for this illustration so that the audio input to the sound 
card will be passed directly to the audio output with a 185 ms delay. Destroy objects and events on exit. 
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Erase RealiIn, ImagIn 


For S = 0 To CAPTURESIZE - 1 Step 2 ‘Copy I to RealIn and Q to Imagin 
RealIn(S \ 2) = inBuffer(S) 
Imagin(S \ 2) = inBuffer(S + 1) 

Next S 


Fig 11—Code for parsing the stereo inBuffer() into in-phase and quadrature signals. This code must be imbedded into the DSP 
subroutine. 


modulation and demodulation of SSB The DirectX runtime driver may be down- You can download this package from the 
signals. Included will be source code loaded from www. microsoft.com/win- ARRL Web www.arrl.org/qexfiles/. Look 
for implementing ultra-high-perfor- dows/directx/downloads/default.asp. for 0902Youngblood.zip. oo 


mance variable band-pass filtering in 
the frequency domain, offset baseband 
IF processing and digital AGC. 


Notes 

1G. Youngblood, AC50G, “A Software- 
Defined Radio for the Masses: Part 1,” 
QEX, July/Aug 2002, pp 13-21. 

2information on both DirectX and Windows 
Multimedia programming can be accessed 
on the Microsoft Developer Network (MSDN) 
Web site at www.msdn. microsoft.com/li- 
brary. To download the DirectX Software 
Development Kit go to msdn.microsoft. 
com/downloads/ and click on “Graphics and 
Multimedia” in the left-hand navigation win- 
dow. Next click on “DirectX” and then 
“DirectX 8.1” (or a later version if available). 
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